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Abstract 



We present the next-to-next-to-leading order QCD analysis of the Gross-Llewellyn Smith 
(GLS) sum rule in deep inelastic lepton-nucleon scattering, taking into account dimension- 
two, twist-four power correction. We discuss in detail the renormalization scheme depen- 
dence of the perturbative QCD approximations, propose a procedure for an approximate 
treatment of the quark mass threshold effects and compare the results of our analysis to the 
recent experimental data of the CCFR collaboration. From this comparison we extract the 
value of the strong coupling constant a™ l (M z ,MS) = 0.115 ± O.OOl(stat) ± 0.005(syst) ± 
0.003(twist) ± 0.0005(sc/ieme). We stress the importance of an accurate measurement of 
the GLS sum rule and in particular of its Q 2 dependence. 
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1. In this paper we continue our investigation of the phenomenological aspects of the 
available next-to-next-to-leading order (NNLO) perturbative QCD approximations to mea- 
surable physical quantities, and in particular of their implications for a precise determination 
of the strong coupling constant a s at different scales (see for a recent detailed discussion 
of this point). In our previous publication we have discussed in detail the NNLO QCD 
predictions for the familiar R-ratios in e + e~ annihilation as well as for the r-lepton decay 
rate. The inclusion of the NNLO corrections calculated in || proved to be quite impor- 
tant, in particular for the latter quantity, as it significantly decreased the sensitivity of QCD 
predictions to the well-known renormalization scheme (RS) ambiguity ||. This makes the 
r-lepton decay rate a suitable place for testing QCD (for more detailed discussions of this 
subject see refs.0,0). Using the results of our work the RS ambiguities were considered 
also in || for the case of the e + e~ annihilation in the resonance region and in |J for the 
r-lepton decay ratio. 

Besides the above-mentioned two quantities, the NNLO calculations are available also 
for some of the deep- inelastic scattering sum rules 0, Efl. The corresponding next-to-leading 
order (NLO) calculations can be found in |)], nj. In this paper we concentrate on the 



Gross-Llewellyn Smith (GLS) sum rule and after discussion of various kinds of theoretical 
ambiguities related to it, present a phenomenological analysis of recent experimental data 
from the CCFR Collaboration |11| . 



2. The quantity of interest in our case is the non-trivial part A of the GLS sum rule 

1 r 1 



GLS = - [ dxF* p+up (x, Q 2 ) 
2 jo 



defined as 

A = (3 - GLS)/3 (2) 

As A = A(Q 2 ) depends, through the Q 2 dependence of the structure function F 3 itself, 
on Q 2 , it would be very useful to have this quantity measured in a broad range of Q 2 
values. Unfortunately the actual situation with the GLS sum rule is more complicated, for 
experimental as well as theoretical reasons, and its analysis burdened with several problems 
related to the treatment of the low x region. The point is that although F$ p , F% p are functions 
of x and Q 2 only, the integral in (|1]) cannot in practice be evaluated from experimental 
data at any finite primary energy. Indeed, the recent experimental analysis of the CCFR 
gives 
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Aex P ((Q 1 ) = 3 GeV 2 ) = 0.167 ± O.OOQ(stat) ± 0.026(systematic) (3) 
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At fixed Q 2 there is always a minimal value of accessible x, determined by the requirement 
that the total hadronic energy W, related to x and Q 2 by 



1 — x 

x 



W 2 = Q 2L_± +m 2 



does not exceed S, the total lepton-proton centre of mass energy: x min = Q 2 /(S + Q 2 — m 2 p ). 
Small values of x thus require small values of the ratio Q 2 /S and if we want to avoid large 
uncertainties due to extrapolation of experimental data to the low x region, we are forced 
to include in the experimental determination of ([!]) a broad range of Q 2 values, down to a 
few GeV 2 . In this region we encounter the following theoretical complications: 

• The quark mass thresholds. The integration over x in (|TJ) means that the corresponding 
values of W span the whole region S > W 2 > m 2 . As all the NNLO calculations have 
been performed for massless quarks, we face the question of the appropriate treatment 
of quark mass thresholds. As for a given Q 2 the quantity determining the effective 
number of massless quarks to be used in evaluating (H) is W, we plot in Fig.l as 
functions of Q 2 the relative contributions w l {Q 2 ) to the GLS sum rule {i = 3,4, 5) of 
the three intervals 

< x < xi, x\ < x < x 2 , x 2 < x < 1 (4) 

Xi = Q 2 + W 2 -m 2 > < = 1 ' 2 ' Wl = Aml W * = 4m * (5) 

in which the effective rif equals 5,4 and 3 respectively {m h = 4.5 GeV , m c = 1.5 GeV 
are the b and c-quark current on-shell masses Q). The curves in Fig.l correspond to 
three widely used sets of distribution functions: Duke-Owens set 1 fl2| , DFLM set 2 
|p~3| and Morfin-Tung (fit S in MS) ||14|| . It is obvious that the fractions w l are only 
weakly sensitive to the specific choice of quark distribution functions. In evaluating 
the fractions w l we have used the parton model formula 

1 r 1 

GLS — - I (u v (x) + d v {x))dx 



2 Jo 

where u v (x) = u(x) — u(x) and d v (x) = d(x) — d(x) are the standard valence quark 
distribution functions. We see that the integral (|l|) is dominated by rif = 3 only for 
very low Q 2 , while for Q 2 = 3 GeV 2 , which is the average value of Q 2 in the data of ||11|| , 



rif = 4 gives the largest contribution, but with significant admixture of nf — 3, 5. The 
sizeable Q 2 dependence of the fractions w l (Q 2 ) comes predominantly from the explicit 
Q 2 -dependence of the boundary value Xi,i = 1,2 and not from the Q 2 -dependence of 
the distribution functions u v (x, Q 2 ), d v (x, Q 2 ). Clearly, the only entirely consistent way 
of incorporating the flavour threshold effects would be to carry out all the calculations 



2 For m c — 1.35 GeV, used in the CCFR analysis, our results change insignificantly. 
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up to the NNLO for massive quarks. As this is practically impossible to do, we shall 
later formulate a parton model based procedure which we believe should reasonably 
approximate the effects of quark mass thresholds. 

Renormalization scheme dependence. For Q 2 in the region of a few GeV 2 the RS am- 
biguities are expected to be phenomenologically important, as was the case in 0. 
Moreover, for the GLS sum rule the RS dependence in the low Q 2 region turns out to 
depend sensitively on the value of n/. 

Contribution of higher twists. Contrary to the case of the r-lepton decay, which is also 
characterized by a rather low value of the natural scale Q 2 , the results [L6| give 



concrete estimates of the dimension-two, twist-four contribution (for recent theoret- 
ical discussions of the dimension-two contributions to the e + e~ annihilation and the 
r-lepton decay R-ratios see ||L7 , Q). The estimates |Tj| suggest that the twist-four 



contribution is quite sizeable and has therefore to be taken into account. 
3. The estimate of the corresponding twist-four contribution to the GLS sum rule, 



performed in [lb| using QCD sum rules formalism, implies for our quantity A 

A tmst ^ = ^ { -^p±, ((O s » = 0.33GeV 2 (6) 

Although the typical QCD sum rules accuracy is about 30%, we prefer to be conservative 
and threfore put a 50% error bar on the above estimate. These errors are consistent with 
the estimate of the higher twist contribution obtained in [|l(J by two other methods, namely 
vector dominance approximation and the non-relativistic quark model. In it was also 
shown that the results of the latter method indicate the problems in the related bag model 
calculations of |TS), which gave negligibly small values of the twist-four contribution. 
For Q 2 = 3 GeV 2 we thus use 

A twtst . 4 = 0.032 ± 0.016 (7) 

Its central value is about 1/3 of the leading perturbative correction in the MS RS.The best 
way to detect the presence of the higher twist contributions would be to study the Q 2 
evolution of the GLS sum rules. Experimentally this is, however, difficult to do and all data 
available so far therefore correspond to the averages over rather broad intervals of Q 2 . 

4. The perturbative part of A(Q 2 ) can be expanded in the renormalized couplant (we 
adopt the notation of |19| ) a(RS) = a s (RS)/7r 



A per ,(Q 2 ) = a(RS)(l + n(Q 2 , RS)a(RS) + r 2 (Q 2 , RS)a 2 (RS) + • ■ ■) (8) 

defined in a particular RS. As indicated in (2) both the couplant a and the coefficients 
do depend on the chosen RS. For the discussion of the RS dependence of physical quantities 
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like (2), the RS may be uniquely denned by the set {a, Ck, k > 2} where c& are related to 
the coefficients of the QCD /5-function. Let us write the couplant a(RS) as a function of the 
renormalization scale variable //, 

da([i, RC) 



din// 



(3(a) = —ba (1 + ca + c 2 a + • ■ ■) 



(9) 



where b = (33 — 2nj)/Q, c = (153 — 19n/)/(66 — 4n/) are the RG invariants while Ci,i > 2 
are free parameters defining the so-called renormalization convention (RC). At the NNLO 
we have two free parameters labelling our RS: c 2 and either a itself or /i, related to it in (|)|). 
The consistency conditions lead to the following explicit relations [HI 



r\ = 61n(/x/A) — p 
r 2 = p 2 - c 2 + (n + c/2) 2 



where 



61n(/i/A) = - + cln 
a 



ca 



\/l + ca + c 2 a 2 1 



+ /(a,c 2 ) 



and 



f(a,c 2 ) 



2c 2 - c 2 



2c 2 



2d 



arctan 



In 



2c 2 a + c 
d 

2c 2 a + c — d 



— arctan • 



2c 2 a + c + d 



-In 



c — d 



d 



4c 2 - c 2 , 4c 2 > c 2 



c + d 



d 



4c 2 , 4c 2 < c 2 



In (^,|11|) p = p(v / Q 2 /A) and p 2 are RG invariants [|1^] and A is defined as 

A = A(2c/b)- c/b 



(10) 
(11) 

(12) 

(13) 
(14) 

(15) 



where A, which is held fixed, is the conventional definition of the QCD A-parameter. Com- 
bining (|I],|r|) we find 



n 



+ cln 



ca 



+ f{a,c 2 ) - p 



(16) 



v Vl + ca + c 2 a? j 

and putting all together we obtain as a function of a,c 2 ,p and p 2 . Note that the energy 

nnl 
pert 



dependence of A""^ enters entirely through the RG invariant p which can be written as 



p = b\n(jQ 2 /A m ) -n(n= JQ 2 , MS) 



(17) 



where we take for the referential RS the MS one. 

The RS dependence of A"^? t (a, c 2 ; p, p 2 ) can therefore be represented by a two-dimensional 
surface in three dimensions. In this picture each point on such a surface represents uniquely 
one RS. Recall that at the NLO Ap l ert was given simply as 



A nl 

pert 



a(2 + ca ln(ca/(l + ca)) — pa) 
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and the corresponding curve was close to a parabola. At the NNLO the surface representing 
A nnl depends non-trivially on the mutual relation of the two RG invariants p and p 2 , and 
in particular on the sign of the latter one. For the e + e~ annihilation and r decay R-ratios, 
P2 < for all rif > 3 and so only this case was discussed in 0]. For ([I]) the situation is more 
complicated as now both cases p 2 < c 2 /4 and p2 > c 2 /4 are physically relevant. 

In H the NNLO calculations of (|l|) were carried out in the MS RS with the result 

r 2 = 41.441 - 8.0272/ + 0.177n 2 (19) 

which coupled with the earlier known formulae for r\ [Kj, and c 2 [50] in the same MS RS 



n = 55/12 -71//3 (20) 

77139 - 15099n f + 325n 2 f , s 

c 2 = f - 21 

9504- 576n/ V ; 

yields an explicit dependence of p 2 (n/) 0- I n particular 

p 2 (3) = 3.438, p 2 (4) = -0.928, p 2 (5) = -5.351 (22) 

To get a quantitative idea of the shape of A"^? t as a function of a and c 2 , for given p and 
p 2 , we look for the stationary points with respect to the variation of a, given by the solutions 
to the equation 

dA nnl 

—E?± = (23) 
da 

In we discussed in detail the situation for p 2 < 0; here we briefly sketch what it looks like 
for p 2 > 0. In this case the solutions of (^) make up a curve in the plane a, c 2 and lie in 
the quadrant a > 0, c 2 > 0. For sufficiently large p, there lies along each such curve a saddle 
point defining the PMS choice of the RS. For large values of the couplant, relevant for the 
IR region, the functional form of the curves, i.e. a(c2, p, P2), can be expressed analytically. 



Expanding (16) in powers of 1/a and keeping only the non-vanishing terms we get, using 



dsn 



Agt = (P2 - c 2 + (7 - P?) a i +( 1 - C --p)a 2 + a + 2( 1 - p) K + 0(l/o) (24) 



where k = l/(3c 2 ) and 



, / /— s 2c 2 - c 2 f n c \ c 

7 = c In(c^) + ^ " -ctan - / === 1 + - (25) 



3 Here we included in P2 the contributions of the light-by- light-type graphs contributing to ri- For the 
phcnomenologically interesting case p > the separation of these contributions in accordance with the 
proposal of [El]] has an entirely negligible effect on the results of this analysis. 
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The equation (p3[) has two simple physical solutions 



P-7 + c/2 + v /3(c 2 -p 2 ) + (p + c/2- 7 ) 2 -3(p-7) 2 

a,i = -^=^=r (2d) 

3(p - 7 + y/c 2 - p2){p - 7 - V c 2 - Pi) 

p - 7 + c/2 - ^(c, - p 2 ) + (p + c/2 - 7 ) 2 ~ 3(p^ 
3(p - 7 + v/c 2 - p 2 )(p - 7 - V c 2 - P2) 
which coincide for c° given as the solution to the equation 

3(c 2 - p 2 ) + (p + c/2 - 7 ) 2 - 3(p - 7 ) 2 = (28) 

At that point 

ai = a 2 = — =^> 00 as p =^> (7 — c/2) (29) 

p - 7 + c/2 

For p < (7 — c/2) only a 2 stays positive and is thus physically acceptable. As we decrease p 
even further the curves move upwards and simultaneously shrink to a point at some c^. The 
lower bound on p follows again from the requirement a 2 > which means 



p > h(c 2 ) = 7(c 2 ) - VC2 - P2 (30) 
The minimum of the function h(c 2 ) lies at oi given as the solution of the equation 

V4c 2 -c 2 (4c 2 - c 2 ) 



- c J4C2 - c 2 + 2c 2 7T - 2 arctan -= - - — n , = (31) 

v V V4c 2 -cV 2Vc 2 - P2 

Substituting the solution of (|3T| ) into (|30| ) we get the lower bound p^f as a function of p 2 . 
For p < Pm^ s , Ap^. is again a monotonous function of the couplant for any c 2 . Solving (|5ID 
for «j = 3 we obtain 



^PMS 



p m '- = 0.863, C2"= lim c 2 PA/5 (p) = 12.29 (32) 

P *Pmin 



It should be stressed that the basic idea of the method of the "effective charges" (EC) 
[Hi ( or the scheme-invariant perturbation theory [p3|| ) should be considered more carefully 
when we go from the NLO to the NNLO. Indeed the condition 

r EC = a EC (33) 

implies at the NNLO merely 

n + r 2 a = (34) 

which has an infinite number of solutions corresponding to the intersection of the surface 
Ap™rt( a ! c 2) with the plane A = a. There are two, one or no intersections for any given c 2 
, depending on the values of p and p 2 . At first glance there is no obvious reason to single 
out one of them. However, requiring r 1; r 2 to vanish separately (as is assumed in ||22|| ) fixes 
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C2 = P2 + (c/2) as well as a. In the following this is what we call the EC result. The equation 
determining the value p^ n at which — ► oo reads 



P™?r, = cln ; = H — — — — f arctan — - — I (35) 

/p 2 + c74 4^ ^ 2^/ 1 J 



and has for nj = 3 the solution 

Pmin(3) = 1-35 

At the NNLO and for n f = 4 the IR fixed points of the PMS and EC NNLO approximants 
are at 

a *PMS = L12; c *pms = _ 2>lgj a *£C = 5>17 C *£C = _ Q 335 ( 3 g) 
while for n/ = 5 

a* PMS = 0.371, cf MS = -10.62, a* EC7 = 0.594 c* 2 EC = -4.95 (37) 

We shall not advocate here any of the popular choices of the RS (PMS, EC or MS) but 
shall take the difference between the PMS/EC and the MS results as a measure of the RS- 
dependence (in the phenomenologically relevant region of p the EC results are practically 
indistinguishable from those of the PMS). As we shall see, typical values of p appropriate 



to the CCFR data |TTJ lie in the region p 6 (2, 6), while p up to 25 might be of interest at 
currently accessible values of Q 2 . 

5. The main difference between the RS dependences for nf = 3,4,5 lies in the IR 
region, which is most of all influenced by the fact that p2 < for nj = 4, 5 and p2 > for 
rif = 3 (see ([H|)) implying no IR stability at the NNLO for rif = 3 even in the PMS and EC 
approaches. The question to what extent is the IR stability, observed for nf = 4, 5 in certain 
schemes, is of physical relevance or merely an artefact of finite order calculations and/or the 
choice of the RS, is difficult to answer on the basis of the available perturbative calculations 
themselves. Nevertheless, if it should really be the case, then a similar behaviour would 
certainly have to be observed for nf = 3 as well. As this is not the case the IR stability 
of the NNLO PMS and EC approximants for nf = 4,5, characterized by fl3"B"D , (|3"TD , has 
probably little physical relevance. Moreover, due to a rather small magnitude of p 2 , A"^f t 
has IR zeros of the PMS and EC approximants at so large values of the couplant that the 
NNLO approximations can hardly be trusted. In the following we shall discuss only the 
region of positive p where the problem of asymptotic explosion does not arise. 

In Fig. 2 we plot the dependence A^tip) on P f° r au three RS: PMS/EC and MS and 
for nf = 3, 4, 5. For PMS/EC results we observe sizeable n/ dependence in particular in the 
small p region. This difference is significant up to p ~ 5. On the other hand the MS results 
are nearly n/ independent!. This somewhat surprising effect is a result of nontrivial partial 
compensation between significant n ^-dependences of the coefficients r^, k = 1,2, as given in 
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( |20| , |T9D , and the couplant a in the MS RS, induced by the n/-dependence of the /3-function 
coefficients bk- 

The same facts, but viewed differently, are presented in Fig. 3 where we plot the RS 
dependence, as measured by the difference between the PMS/EC and MS results, for n/ = 
3,4,5. The PMS and EC approaches are practically indistinguishable in our plots. At the 
NNLO and for fixed p the RS dependence diminishes as nj grows, while at the NLO it 
remains about the same. The importance of the NNLO corrections with respect to the NLO 
ones can be assessed from Fig. 4, where we plot the comparison between NLO and NNLO as 
a function of p for all three values if n; and for both PMS/EC and MS RS. We see that 

• in the MS RS we find A^ rt (p) < A™' t (p) for all values of p, the difference A^ t - A^ rt 
being a decreasing function of nf, 

• in the PMS/EC approaches the difference A™g rt — A"™^ is positive for n/ = 5, close 
to zero for nf = 4 and negative for rif = 3. Moreover this difference is sizeable and 
thus of phenomenological relevance. For the consequences of these two observations 
see further discussion; 

• in the phenomenologically important region of Q 2 the incorporation of the NNLO 
corrections into the analysis decreases the difference between the PMS/EC and MS 
results. 

The curves shown in Figs. 2,3,4 can be used for straightforward determination of the QCD 
parameter Aj^g from experimental data on A. For easy use we have fitted them by an 
analytical expression of the form 

AgUp) = £ % (38) 

3=0 P 

in the wide interval p G (2,26), with i labelling one of the 18 combinations of the order 
(NLO or NNLO), RS (PMS, EC or MS) and n f (3,4,5). The values of all the parameters rf 
are given in Table 1. The extraction of A^jg from a given experimental value A exp proceeds 
in two steps: first one solves the equation A p %. t {p) = A exp and then uses the solution pW p in 
the formula 



A Ss = ^ ex P 



V 1 (MS)+pW p ) 



f) (39) 



The analytical parametrization of the dependences A p l l rt (p) is useful also for the deter- 
mination of the error a\ of the extracted Af^g from the error o"a on A exp : 

(40) 



Ams b\dA(p)/dp\ 
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The relation ( |39"D holds for massless quarks and thus makes good sense only provided rif 
is well defined and fixed. As in practice this is usually not the case, we have developed a 
procedure which takes into account at least in an approximate way the quark mass thresholds. 
For all four combinations of the order (NLO and NNLO) and RS (EC/PMS and MS, EC 
and PMS being indistinguishable for our purposes), it consists of the following steps: 

1. From the information on valence distribution functions we determine, for a given Q 2 , 
the fractions w l (Q 2 ) defining the relative importance of the contribution of rif = 3, 4, 5 
to ([l]). We do this by taking the average of results corresponding to the three mentioned 
sets of quark distribution functions. For Q 2 = 3 GeV 2 we find: w 1 = 0.20, w 2 = 
0.51, w 3 = 0.29. 

2. From (j39"D A^ appropriate for rif massless quarks is determined. 



3. Using the formulae from [24j these values are then translated to A^-(rif) where the 



argument n f keeps track of the source of this resulting A- 



MS' 



4. The average value A^ = Yn=a u, *-^o5 , (0 similarly for errors) is evaluated. 

5. Using the average value A^ we go back and calculate the corresponding values of 
A (n/)_ 

MS 

6. From these values of A^j- we calculate, according to (|17D , the corresponding p{rif) and 
then finally evaluate the averages A pert = £)f =3 w l A pert {pi). 

For this procedure to be self-consistent we should come at the end of step 6 close to 
the experimental value of A exp |3|. If the data do not correspond to a fixed value of Q 2 the 
whole procedure should be folded with the known Q 2 dependence of the data. As this is 
difficult to do we have carried out the above procedure for fixed Q 2 = 3 GeV 2 , equal to 
the average value of Q 2 in the data of [1TJ. We have included the higher twist contribution 



as given in (|7|) Q but carried out the whole analysis also for the case of no higher twists at 
all. In Table 2 we present in detail the results of our procedure for the case of the central 
value of Atwist-4 = 0.032. Notice that the averages in the last column are indeed quite 
close to A exp — 0.032 = 0.135. Similar tables have been constructed for the upper (0.048) 
as well as lower (0.016) estimates of twist-4 contributions and also for the case of no twist-4 



Table 3. From Tables 2 and 3 we conclude that: 



contribution at all. In these latter cases only the extracted values of A^ are presented in 



4 We neglected in our analysis the kinematical power corrections and the twist-six corrections since it was 
shown in ref. that they are significantly smaller than the central value of the twist-four contributions (0) . 
It might be interesting to take these effects into account in order to compare their value with the assumed 
by us 50% error bars of the twist-four contibutions. 
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• the inclusion of the higher twists is quite important as it lowers the central value of 
A m by about 100 MeV; 

• in the MS as well as PMS/EC approaches A^(NNLO,n / = 3) < A^(NLO,n/ = 3) 
while in the PMS/EC approach they are almost the same (smaller) for rif — 4 (n/ = 5). 
This is a consequence of the first two observations of Sec. 5, namely that Ap l ert < A""/ t 
in the MS RS and for PMS/EC with n/=3, while in the latter case for rif =4 (rif— 5) 
we have Ag rt « A™', (A^ ri > A^) (see Fig.4); 

• the inclusion of the NNLO corrections reduces substantially the RS dependence as 
measured by the difference between the results for PMS/EC and MS approaches while 
its influence on A^jg is smaller than the effects of the twist-four corrections; due to 

• the errors of the extracted A^ values are dominated on one side by the systematical 
experimental uncertainties and on the other by uncertainty in the higher twist contri- 
bution. The reduction of the systematical experimental errors (which are the biggest 
single source of errors and make about half of the total error), to the level of statis- 
tical ones would make the GLS sum rules a very good place for testing perturbative 
QCD provided a more accurate estimate of the higher twist were to be available. But 
even taking into account the combined effect of all discussed errors the accuracy of the 
extracted A^ is not bad; 

• the results of the averaging procedure are not far from those which correspond to fixed 
rif = 4. This is due to the fact that for Q 2 = 3 GeV 2 rif = 4 is dominant while 
n/ = 3, 5 contribute with comparable strength, thus roughly balancing each other. For 
other Q 2 the situation may be quite different. 

6. Let us now summarize our results. At the NLO and in the MS RS we find 

Agl = 317 ± 23(stat) ± 99(syst) ± 62(twist) MeV (41) 

while the NLO PMC/EC result reads 

A^ = 241 ± U(stat) ± 60{syst) ± Activist) MeV (42) 

where the third error, due to higher twist contributions corresponds to the limits given in 
fl7|). Notice that at the NLO the value of A^ extracted within the PMS/EC approach is 
smaller than the one obtained in the MS analysis. A similar situation occurs in the NLO 
analysis of the BCDMS data on F 2 (x,Q 2 ) without higher twist contributions [^] and in 
the NLO analysis of the BCDMS and SLAC data with higher twist contributions [21]]. Our 



analysis differs from that of [26| and other similar ones in the way higher twists are treated. 



While we use concrete estimates for twist-four contributions (though with sizeable errors), 
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previously they were simply parametrized by the corresponding free parameters which were 
fitted together with pure perturbative expressions from experimental data. 

At the NNLO the difference between the MS and the PMS/EC results becomes consid- 
erably smaller: 

MS: A^| = 265±l8(stat)±80(syst)±50(twist) MeV (43) 
PMS/EC : A^ = 249 ± IQ(stat) ± 70(syst) ± Ah(twist) MeV (44) 



For a careful reader we add the following comment concerning the different relation between 
the errors at the NLO and the NNLO in PMS/EC and MS RS. As follows from Figs. 2,3,4, 
in the PMS/EC approaches the central NNLO value of A^ is slightly larger than the NLO 
one. Furthermore, for nf=4 the derivative of the PMS/EC approximant to A(p) is somewhat 
bigger at the NLO than at the NNLO. This leads (see (|4T)|)) to slightly lower error <j\ at the 
NLO. The combination of these two small effects results in a slight increase of both theoretical 
and experimental uncertainties of A^ at the NNLO in PMS/EC approach compared to the 
corresponding NLO analysis (compare (43) with (41)). 
Combining the MS and PMS/EC results we finally get 



A^|(NLO ) = 279 ± 19{stat) ± 80(syst) ± 50{twist) ± 38(scheme) MeV (45) 
A^(NNLO) = 257 ± 17(stat) ± 75(syst) ± A7{twist) ± 8(scheme) MeV (46) 



Notice that with respect to the NLO results the NNLO statistical, systematical and higher 
twist uncertainties are slightly smaller, while the error due to the RS-dependence is reduced 
significantly. This is the main effect of including the NNLO corrections in the analysis of 
GLS sum rule. Our results are in very good agreement with those of the detailed NLO 
analysis |27j (carried out in the MS RS and including the phenomenological parametrization 
of the higher-twists effects) f] of combined BCDMS and SLAC data on F2(x,Q 2 ) structure 
function: 

A^| = 263 ± 42(exp) MeV (47) 

As, however, our analysis suggests the importance of including the NNLO corrections, it 
would be very desirable to investigate the influence of the yet uncalculated NNLO corrections 
to F 2 (x, Q 2 ) on the result ( |47| ) as well. 

So far our analysis has been aimed primarily at the determination of A^ but it tells 
us also something about the role of the higher twist contributions to the GLS sum rule. 
Indeed, Table 1 suggests that if the twist-four contribution (|7|) is ignored or its magnitude is 



negligibly small (like in the bag model calculations |T8[), the agreement between the results 



of our analysis and those of |27[ worserns. This lends a posteriori support to the QCD sum 



5 The recent NLO MS and scheme-invariant fits [§(J of the BCDMS and SLAC data, including the higher 
twists, gave smaller values of Aj^g. It would be interesting to understand the origin of these smaller values. 
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rules estimate [IB| of the higher twist effects. An independent and more accurate estimate 
of these higher twists, using for instance lattice methods (for a review of such calculations 
see |[28||), would, however, be very welcome. This analysis might also shed the light on the 
role of other possible 1 / Q 2 -contributions, recently discussed, from a purely theoretical point 
of view, in [llTL pi. We intend to consider related problems in the future. 



The above results on imply for a s at the scale y/Cp = \/3 GeV the following values 



MS 

, nl 



af(^jQ 2 , MS) = 0.318 ± OMO(stat) ± 0.0A2(syst) ± 0.026(twist) ± OMO(scheme) (48) 

as nl {\/0*, MS) = 0.315 ± O.OlO(stat) ± O.OU(syst) ± Q.Q28(twist) ± 0.005(sc/ieme) (49) 

where in accordance with the discussions of Sec. 2 we have taken rif = 4. The errors a as of 
the above results on a s are related to the errors a a of (44), (45) as follows 

where F is the r.h.s. of eq.(10). 

Finally, to facilitate easy comparison with other determinations of a s we have evolved it, 



using the fomulae of Pl |, and eqs.(^8|),(^9|), through the b-quark threshold up to M z and 
obtained 

a"'(M z ,MS) = 0.116 ± O.OOl(stat) ± 0.005(syst) ± 0.003(twist) ± 0.002(sc/ieme) (51) 

a™ l (M z , MS) = 0.115 ± O.OOl(stat) ± 0.005(syst) ± 0.003(twist) ± 0.0005(sc/ieme) (52) 
These results can be compared with 

a s (M z ,MS) = 0.113 ±0.003(ea;p) ±0.004(t/ieor) (53) 

obtained in the NLO analysis in the MS RS P7[] , where the standard two-loop inverse-log 
approximation for a s was used. The theoretical error comes from the NLO estimates of the 
RS and factorization-scheme uncertainties []. 

Our results also agree with the recent world average |I[ 

a s (M z ,MS) = 0.112 ±0.007 

determined on the basis of the results of various deep-inelastic scattering experiments. The 
quoted error bars correspond to the combination of experimental and theoretical uncertain- 
ties with the dominant role of the latter ones. This average has in turn been used in the 
determination of the overall world average |l| 

a s (M z ,MS) = 0.118 ±0.007 



6 We do not use this approximation but are solving ( jl2| ) numerically to get a s as a function of the scale 
M = \fO^- Within this procedure, we got at the NLO 0.115 instead of the above mentioned result 0.113. 
This theoretical difference lies within the total theoretical error of the result (|53|). 
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with the typical error (mainly theoretical one). This world average sums up all available 
results on a s , including the relatively larger values coming from the analyses of the R-ratio 
in the e + e~ annihilations into hadrons and of the total width of Z° decay into hadrons, both 
of which, however, have rather large experimental errors (see |29|, |30], [L|). 

7. In this work we have presented the NNLO QCD analysis of the GLS sum rules 



using the new data of the CCFR Collaboration |]TTJ . The inclusion of the NNLO corrections 
substantially improves the situation as far as the RS dependence is concerned and makes 
this quantity a potentially good place for testing the perturbative QCD. Our result 

a™ l {M z ,MS) = 0.115 ± O.OOl(stat) ± 0.005(syst) ± 0.003(twist) ± 0.0005(sc/ieme) 

is in good agreement with the current world average. However, three problems, two theoreti- 
cal and the other experimental, stand, in the way of a more precise determination of a s [Mz) 
or Ajyfg from this data. The first one concerns the treatment of flavour thresholds and we 
have suggested an approximate procedure to take this effect into account while still using 
the calculations for massless quarks. The second problem is connected with the necessity of 
getting more precise estimates of the theoretical errors of the higher twist terms. The third 
problem is related to the experimental errors, which are dominated by the systematical ones. 
These errors are currently larger then the theoretical errors of the higher twist terms and 
any attempt to reduce them would thus be very usefull. It would also be very interesting to 
have accurate data on (|1]) as a function of Q 2 in a broad range of Q 2 values. Their analysis 
would certainly lead to better understanding of the role of the higher twist contributions to 
deep inelastic scattering at low and moderate Q 2 . 

We intend to carry out an analysis similar to the one presented in this paper also for the 
Bjorken sum rule for polarized electroproduction structure functions, since the appropriate 



NNLO perturbative calculations M as well as power correction estimates [31], inspired by 



the analysis of f32| , are available and the relevant experimental data will soon be obtained 



by CERN and SLAC groups. 
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Figure captions 



Fig.l: The fractions w\i = 1,2,3 as the functions of Q 2 for three widely used sets of quark 
distributions functions. 

Fig. 2: The dependence of A(p) from n/ for the PMS/EC approaches and the MS RS (denoted 
in this and allthe following figures by MSb), at the NLO and the NNLO. 

Fig.3: The RS dependence of A(p) for n f = 3, 4, 5 at the NLO and the NNLO. 

Fig. 4: The comparison between the NLO and the NNLO results for A(p) in the cases of 
71/ = 3, 4, 5. 

Note: In the process of the study of Fig.l-Fig.4 the readers should cnange MS to MS every- 
where. This was not done in view of the technical problems in dealing with the corresponding 
PS-file. 
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n f 


(i) 


(i) 
r l 


(i) 
r 2 


(i) 


(i) 

r 4 


(i) 
r 5 






3 


0.008137 


0.645500 


0.133457 


-0.413737 


-3.100763 


7.573763 




PMS 


4 


0.008448 


0.653533 


0.076082 


-1.053512 


-3.038260 


7.176102 






5 


0.004056 


0.814661 


-1.174005 


0.267708 


2.229508 


-2.753092 


N 




3 


0.007837 


0.653591 


0.081509 


-0.416852 


-2.659066 


7.374106 


N 


ECH 


4 


0.007270 


0.691163 


-0.264749 


-0.428416 


-1.482951 


3.615289 


L 




5 


0.003300 


0.834586 


-1.291855 


0.353730 


2.839872 


-3.593072 






3 


0.004927 


0.760169 


-0.963719 


0.123290 


2.019464 


-2.445137 




MS 


4 


0.008325 


0.660719 


-0.062042 


-1.260045 


-2.644704 


7.022455 






5 


0.008725 


0.667920 


-0.016647 


-1.436806 


-2.770146 


7.534455 






3 


0.004645 


0.761183 


-0.904944 


0.435364 


2.299329 


-3.332708 




PMS 


4 


0.008607 


0.648385 


0.155498 


-1.095457 


-3.068545 


7.277303 






5 


0.006449 


0.732606 


-0.245552 


-0.864299 


0.201203 


2.163065 






3 


0.004908 


0.754859 


-0.903029 


0.442944 


1.925275 


-2.829597 


N 


ECH 


4 


0.004350 


0.782570 


-0.926641 


0.446002 


2.117838 


-3.104819 


L 




5 


0.008235 


0.675962 


0.148810 


-1.066902 


-3.437890 


7.772996 






3 


0.005703 


0.745372 


-1.420683 


0.564036 


3.195462 


-4.264516 




MS 


4 


0.002404 


0.851741 


-2.019862 


2.597744 


-0.376554 


-1.574096 






5 


0.008854 


0.668406 


-0.247080 


-1.211435 


-2.402028 


6.958375 



Table 1: The values of the coefficients for the i — th option where the option is defined 
by the combination of the order, RS and rif. 
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T 


MS 
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4 35 


9^0 


267+1 8(79") 
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140 


1 37 






5 


4.45 


317 




0.121 








3 


4.50 


212 




0.147 






PMS 


4 


4.70 


230 


240±14(60) 


0.140 


0.136 


N 




5 


5.00 


278 




0.122 




L 




3 


3.10 


300 




0.139 






MS 


4 


3.50 


307 


318 ±23(99) 


0.138 


0.135 






5 


4.00 


352 




0.127 





Table 2: Results of the procedure described in the text for the case with twist-four contri- 
bution given by the central value of the estimate (5). The values of the A-parameter are in 
MeV and the first (second) error corresponds to the statistical (systematical) experimental 
errors. 





RS 


HT=0 


HT=0.016 


HT=0.032 


HT=0.048 


NNL 


PMS 


333 ± 15(65) 


294 ± 15(65) 


250 ± 16(67) 


204 ± 16(70) 


MS 


360 ± 19(83) 


313 ± 19(81) 


267 ± 18(79) 


215 ± 17(74) 


NL 


PMS 


319 ± 13(55) 


281 ± 13(55) 


240 ± 14(60) 


201 ± 14(60) 


MS 


435 ± 20(87) 


378 ±21(90) 


318 ±23(99) 


254 ±24(105) 



Table 3: The values of A^- (in MeV) as extracted from data of the CCFR Collaboration 
under different assumptions of the higher twist (HT) contribution. The errors are the same 
as in the Table 2. 



18 



